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A class of spherical functions is studied which can be viewed as the matrix generalization of Bessel 
functions. We derive a recursive structure for these functions. We show that they are only special 
cases of more general radial functions which also have a, properly generalized, recursive structure. 
Some explicit results are worked out. 



I. INTRODUCTION 



O ' In 1957, Harish-ChandraE] derived a famous formula for a certain class of group integrals. Let Q be a compact 
' semi-simple Lie group and let a and b be elements of its Cartan subalgebra H, then 

■ f , / i rTi\ 1 exp (tr u>(a)6) , „. 

' W)exp(trt/-W&) = — V -g— > . (1.1) 



Jues " v ' " v ' |W| ^ n(a)nM&)) 

! Here, dfi(U) is the invariant measure, 11(a) is the product of all positive roots of TL and W is the Weyl reflection 
■ group of Q with \W\ elements w. 

This result depends crucially on the condition that a and b are in the Cartan subalgebra TC. In other words, 
U~ 1 aUb has to be in Q . If one re place s a and b in the integral on the left hand side with more general matrices x 



and k which are not in 7i, formula (1.1) is not valid anymore. The spherical functions introduced by Gelf andBEj form 



an important class of sach integrals which are, in general, not covered by Harish-Chandra's result ( |l . l[) . In another 
work, Harish-ChandraE2l studies in great detail the harmonic analysis involving these spherical functions. In a more 
physics oriented contribution, Olshanetsky and PcrelomovEil discussed them in the framework of quantum integrable 
• systems. 

Here, we wish to address spherical functions of the following kind: we take x and k as diagonal matrices containing 
the eigenvalues of a Hamiltonian in a matrix representation. The Hamiltonian is diagonalized by the integration matrix 
U. In particular, we assume that the Hamiltonian U~ 1 xU or, equivalently, UkU^ 1 is real-symmetric, Hermitean or 
Hermitean self-dual. Thus, Q is the orthogonal, the unitary or the unitary-symplectic group. We will refer to these 
spherical functions as matrix Bessel functions. We notice that the uni tary case is special: since it so happens that the 



'"^ , eigenvalues x and k do lie in the Cartan subalgebra 7i, the result (LI) applies and coin cides with the Itzykson-Zuber 
d • formulate. In the orthogonal and the unitary-symplectic cases, however, formula (LI) is not valid. 

We choose the term matrix Bessel function for the spherical functions to be discussed here, because they can 
be viewed as a natural extension of the ordinary vector Bessel functions. However, due to the rich features of the 
spherical functions, other extensions relating to ordinary Bessel functions are equally natural. Related functions 
have beettL. discussed aniL terms similar to matrix Bessel functions have akpady been used by HertzES, by Gross and 
KunzcE3Ej, by Holmaro ancl-by Okounkov and 01shansk£3. KontsevictO introduced the matrix Airy functions. 

Duistermaat and Heckmanu developed- a stationary phase approach involving localization for a class of spherical 
functions, see also the treatise by SzaboEi 

Remarkably, our matrix Bessel functions are only special cases of more general objects which we call radial functions. 
Moreover, there is an important connection to the Calogero-Sutherland models which we will discuss separately, see 
below. 

The matupL. Bessel functions are of considerable interest for applications in physics. They appear in Random Matrix 
TheorycJ'lSEj which models spectral fluctuations mL complex systems, such as quantum chaotic ones. In particular, 
they are the kernels of Dyson's Brownian motionalZI describing crossover transitions between different symmetry or 
invariance classes. Unfortunately, only the case of broken time-reversal invariance can be treated explicitly with the 
help of the Itzykson-Zuber formula. In the physically important cases of conserved time-reversal invariance, the 
kernels are not known analytically, as argued above. Muirheado discusses spherical functions in the framework of 
multivariate statistical theory. In his book, an expansion in terms of Jack polynomials for the orthogonal case can be 
found. Such an expansion for arbitrary Dyson index was recently worked out by Okounkov and OlshanskiEl 
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The goal of the present paper is to explore the structure of the radial functions which contain the matrix Bessel 
functions as special cases. In particular, we show how explicit results can be obtained. The paper is organized as 
follows. In Sec. |Tl], we briefly review some properties of the vector Bessel functions. In doing so we wish to help the 



reader in developing an intuition for the matrix Bessel functions which wc introduce in Sec. [II. In Sec. IV, we state 
and derive a fundamental recursive structure for matrix Bessel functions. We show in Sec. |\^that this recursion is 
an iterative solution of general radial functions which contain group integrals defining the matrix Bessel functions 
as special case. Sees. IV 
and explicit formulae 
separately in Sec. VII. In Sec 
the appendix. 



and are our main results. In Sec. [VD we illustrate how the recursion can lead to closed 
Because of its special importance, we discuss the connection to Calogero-Sutherland models 
we summarize and conclude. Various aspects and calculations are collected in 



II. VECTOR BESSEL FUNCTIONS REVISITED 

Before turning to the matrix case, we compile, for the convenience of the reader, some well known results for the 
vector case. 

In a real, d dimensional space with d = 2, 3, 4, . . ., we consider a position vector r — (x\, . . . , Xd) and a wave vector 
k = (fci, . . . , kd). The plane wave exp {ik ■ satisfies the wave equation 

Aexp {ik ■ fj = — k 2 exp {ik ■ r\ (2-1) 
where we define the Laplacean as in the physics literature, 

Q2 d Q2 

A = ^ = £&r (2 - 2) 

i=l 1 

The zeroth order Bessel function in this space is the angular average of the plane wave, 

X {d \kr) = J dfl exp {ik ■ rj , (2.3) 

over the solid angle Q, defining the orientation of either for fc. In our context, it is advantageous to take Q as the solid 
angle of k. Obviously, only the relative angle between r and k matters and x {kr) can only depend on the product 
of the lengths r = \ f] and k — \k\ of the two vectors. We normalize the measure dn with the volume 2n d / 2 /T(d/2) of 
the unit sphere, i.e. we have 

J dQ = 1 . (2.4) 

Thus, by construction, we also have 

X {d) (0) = 1 . (2.5) 

It is convenient to view r as the azimuthal direction of the coordinate system in which we measure ft. Thus, in these 
spherical coordinates, one finds k ■ r — fcr cosi? where § is the azimuthal angle. The measure dfi contains sin d ~ 2 $ 
and one has 



yVr((d- l)/2) J 



^- 2)/2 W2) 'f^r (2-6) 



where J v {z) is the standard Bessel functionS of order v. The functions (2.6) are often referred to as zonal functions. 

There is a remarkable difference for the functions X (^ r ) if one compares even and odd dimensions. For example, 
one has in d = 2 dimensions X {kr) = Jo(kr) and in d = 3 dimensions x^ 3 ^(A:r) = [t: /2) 1 / 2 Ji/ 2 {kr) / (kr) 1 / 2 = jo(kr) 
with the spherical Bessel function jo(z) of zeroth ordera. In d — 2 dimensions, Jq(z) is a complicated infinite series 
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in the argument z, in d = 3 dimensions, however, jo(z) is the simple r atio jo(z) = sin z/z. One easily sees how this 
generalizes. Upon introducing £ = cost? as integration variable in Eq. (2J3), one finds the representation 



ViT((d-l)/2) 
In dimensions d > 3, this can be cast into the form 



£\xp(ikr0(l-ef- 3)/2 d!;. (2.7) 



y W (hr) = 2F ( rf / 2 ) V ( (d " 3)/2 ^ ^ Sinfcr (2 8) 

X K ' y^T((d- l)/2) ^ V A* J d(kr) 2 » kr ' V j 

For eve n ri, the exponent (d— 3)/2 is a fraction —1/2, +1/2, +3/2, . . ., and the function (l — £ 2 )' d 3 ^ 2 in the integrand 
in Eq. (2/7) is an infinite power series. This yields, for d = 4, 6, 8, . . ., the complicated power series (2^) involving 
an infinite number of inverse powers of kr. However, if d is odd, the exponent (d — 3)/2 is an integer 0, 1,2,.. ., and 

the function (l — £ 2 ) 3 ^ 2 is a finite polynomial of order (d — 3)/2 in £ 2 . Thus, x^ d \kr) acquires a comparatively 
simple structure, because it only contains a finite number of inverse powers of kr. Formally, this means that for odd 
d all binomial coefficients for [i > (d — 3)/2 are zero. 

The differential equation for the functions x^ d \kr) is easily obtained by averaging Eq. (2.1) over the solid angle fl 
of k, i.e. by integrating both sides, 

A J dfiexp(ifc-r) = -k 2 J dnexp(ik-i^j . (2.9) 

We notice that the Laplacean A commutes with the integral, because the former is in the space of the position vector, 
the latter in the space of the wave vector. Moreover, the integral trivially commutes with k 2 = k 2 . Hence, one arrives 
at 

A rX {d) {kr) = -k 2 X {d) {kr) . (2.10) 
Since x^ d \kr) depends exclusively on radial variables, we replaced the full Laplacean A with its radial part 

1 d rf i d d 2 d - 1 d , n , 



1 1 dr dr dr 2 r dr 

Ad) 



In general, there are two fundamental solutions x+ (kr) and x~ (kr) of the differential equation ( 2.1C ) which behave 
as exp (±ikr) / (kr)^ d 1 '^ 2 for large arguments kr. Thus, to obtain the full solutions, one can make the Hankel ansatz 

X^ikr) = e -^^^\kr). (2.12) 



(kr) 



Here , w+ \kr) is a function with the property w±\kr) — > 1 for kr — > oo. The differential equation follows easily from 
Eq. (2.11) and is given by 

d 2 , d d- \ (d-\ \ 1 \ {d) 



For d > 3, one uses the ansatz as an asymptotic power series 



which yields a recursion for the coefficients 

1 / , , d-\ (d-1 



i2(n+ 1) 



K» + 1) s- -5- ~ 1 «M . ( 2 - 15 ) 



with the starting value ao = 1 . A special situation occurs when the integer running index fi reaches the critical value 
fi c = (d — 3)/2. If d is odd, /i c is integer and the recursion terminates at [i = /x c , i.e. one has a M = 0, fi > fi c . Thus, 
the asymptotic series becomes a finite polynomial in inverse powers of kr. However, if d is even, /j, c is half-odd integer 
and the series cannot terminate, it is always infinite. This explains the different structure of the Bessel functions in 
even and odd dimensional spaces from the viewpoint of the differential equation. 

In App. ^ we discuss an alternative integral representation which has an interesting analogue in the matrix space. 
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III. MATRIX BESSEL FUNCTIONS 



We compile the basics features of the m atrix s paces we want to work with in Sec. [II A, before we define the matrix 



Bessel functions as group integrals in Sec. [II B 



Two general aspects are shifted into the appendix. First, we present an interesting alternative integral representation 
in App. [b]. Second, the matrix Bessel functions play a crucial role in harmonic analysis or, equivalently, in Fourier= 
Bessel analysis in matrix spaces. For the general theory, we refer the reader to Harish-Chandra's treatise in Ref.cJ 
and to Helgason's bookcil. However, to achieve our goal of being explicit, we collect, for the convenience of the reader, 
some results for the Fourier-Bessel analysis of invariant functions in matrix spaces in App. |c[ 



A. Basics and Notation 



We introduce N x N matrices H whose elements H nm , n,m = 1, . . . , N are real, complex or quaternion variables. 
In other words, each element H nm has (3 real components H^l, a = 0, ...,(/?— 1) with (3 = 1, 2, 4, respectively, 

H nm = J2 H n2r {a) . (3.1) 



Here, we use the basis T^ a \ a = 0, . . . , (/3 — 1). We have — 1 for the real case with (3 = 1. For the complex case 
with j3 = 2, we have = 1 and = i. Finally, we have 



,(2) _ 



1 
1 

-i 
-i 



-(i) _ 



-O) _ 



+1 
-1 

' +i 
-i 



(3.2) 



in the quaternion case for (3 = 4 where the r^ a \ a = 1,2, 3 are the Pauli matrices. We notice that the total H is 
a 2N x 2N matrix for (3 — A. However, here and in the following, the dimensions that we use always refer to the 
number of matrix elements such as H nm . These are scalar for (3 = 1,2 and quaternion for (3 — 4. The label (3 is often 
referred to as Dyson index. 

We assume that the matrix H is real symmetric, Hermitean or Hermitean self-dual in the three cases (3 = 1, 2, 4. We 
always write ift = H to indicate this symmetry. There are N independent real variables H nn = ff^n , n= 1, . . . , N 
on the diagonal and f3N(N — l)/2 independent real variables ffnm, a = 0, ...,(/? — 1), 1 < n < m < N outside the 
diagonal. We write the volume element of H in the form 



N 



0-1 



(3.3) 



71=1 



The matrix H is diagonalized by the matrix U, with columns U n , n = 1, . . . , Un- Depending on the value of /3, the 
matrix U is either orthogonal, unitary or unitary-symplectic. Following Gilmore's notationEil, we write U € U (N; (3) 
with U(N; 1) = SO(N), U(N; 2) = U(N) and U{N; 4) = USp(2N). The volume of these groups is given by 



volU(N;f3) 



We use it to normalize the invariant measure dn(U) of U £ U(N; (3) to unity, 

dn{U) = 1 . 



(3.4) 



(3.5) 



The N real eigenvalues x n , n — 1, . . . ,N of H are ordered in the diagonal matrix x. We have x = diag (xi, . . . , xn) 
for /3 = 1 and (3 = 2. For (3 = 4, the eigenvalues are doubly degenerate and we have x = diag (xi, xi, . . . , xn, Xn)- 
Physically, this doubling of the eigenvalues is due to Kramer's degeneracies. Thus, the diagonalization reads 



H = U^xU , 



with 



U^xUryi 



(3.6) 
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The, diagonalizing matrix U has the property U 1 = W . The volume element in eigenvalue-angle coordinates is given 
by§H 

d[H] = C ( f\A N (x)fd[x]dfx(U) (3.7) 
where d[x] denotes the product of all differentials dx n . We have introduced the Vandermonde determinant 

Ajv(s) = [l^n-Sm) ■ (3.8) 



n < rn 



The normalization constant 



7r /3AT(A r -l)/4 r JV (j g/ 2 ) 



obtains from the constants given in Mchta's bookEZI and from Eq. (3.4). 

To avoid inconveniences and to ensure a compact notation, we define the trace Tr and the determinant Det with 
Tr = tr and Det = det for j3 = 1, 2 and with 

1 



TrK = -tiK and Det If = VdctK (3.10) 

in the case /3 = 4 for a matrix K with quaternion entries. If k denotes the diagonal matrix of the eigenvalues of a 
real-symmetric, Hermitean or Hermitean self-dual matrix, it is also useful to define the associate matrix k. In all 
three cases (3, it is the N x N matrix k = diag (fci, . . . , /cat), i.e. we have k = k for (3 = 1, 2 and no degeneracies for 
= 4. 



B. Integral Definition and Differential Equation 

As in the case of vector Bessel functions, we start in the matrix case with the plane wave. For two matrices H and 
K with the same symmetries = H and — K, we introduce the matrix plane wave as exp (iTr H K) where the 
trace is the proper scalar product in the matrix space. The matrix plane wave has the property 

{2v) N v mN-i),2 J d[H]ex V (mHK) = 8{K) (3.11) 

where 8(K) is the product of the 5 distributions of all independent variables. We define the matrix gradient d/dH 
and the Laplacean operator 



A = Tr 



E3 + oEE77^ ( 3 - 12 ) 



QJJ2 p)TT<,0)2 2 ^ ^ f)n^ 2 

71=1 un nn n<ma=0 urL nm 

which acts on the matrix plane wave as 

A exp (zTr HK) = — Tr K 2 exp (iTr HK) . (3.13) 

We notice that, for (3 = 4, inconvenient factors of two would occur if we used tr instead of Tr . 

Analogously to vector Bessel functions, we define the matrix Bessel functions as the angular average 

^\x,k) = J dfi(U)exp(iTvHK) . (3.14) 



The diagonal matrix k contains the eigenvalues of K which is diagonalized by a matrix V such that K = V'kV 

N 



Due to the invariance of the measure dfi(U), the matrix V is absorbed and the functions <S>x(x,k) depend on the 



eigenvalues x and k only, 

§%\x,k) = I dLi(U)ex-p(iTrU' i xUk) . (3.15) 
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Thus, in the scalar product Tr HK , solely the relative angles between H and K matter. The matrix Bessel functions 
are symmetric in the arguments, 



(3.16) 



and normalized to unity, 



*$\x,0) = 1 



and 



$^(0,fc) = 1 



(3.17) 



due to Eq. (3.5). These are spherical functions in the sense of Ref.t 



As in the vector case, the differential equation is obtained by averaging Eq. ( 3.13 ) over the relative angles, 

A J dfx(V)exp(iTrHK) = -Tr K 2 J dfx(V) exp (iTr HK) . (3.18) 

Again, the Laplacean A commutes with the integral, because the former is in the space of the matrix H , the latter 
over the diagonalizing matrix V of K. The integral also commutes with Tr K 2 = Tr k 2 . Due to the symmetry between 
H and K, the integral is obviously identical to the definition (3.15) and we find 



A x ^\x,k) = -Trk?9%>(x,k) 



(3.19) 



Since the matrix Bessel function <$>(£\x, k) depends only on the radial variable s, i. e. on the eigenvalues, we replaced 



the full Laplacean with its radial part A x . Because of the transformation rule (3.7), it reads 



N 

E 



1 9 p d 

< IAjvW^; 1 N[X)l dx~n 



TV 



& 



E Qx 2 + E 







d 

dx r , 



dx n 



(3.20) 



We notice that these steps are fully parallel to the corresponding discussion in Sec. [□]. Importantly, due to the 
symmetry ( 3.16| ), the functions <f>x\x, k) must also solve the differential equation in the k n , n — l,...,N which 
results from Eq. (3.19) by exchanging x and k. Obviously, this a very restrictive requirement. 

Comparing the radial operator (2.1C) in the vector case and the radial operator ( |3.20] ), we see that it is the j3 that 
corresponds to the spatial dimension d or, more precisely, to d — 1. The role played by the matrix dimension N is 
a different one. To illustrate this, we study the two simplest cases. First, we can formally set N = 1 and find from 



the definition (3.14) that <fr\ (x,k) = exp(ix\ki) where H\\ = x\ and K\\ = k\. In this case, the matrix U has 
dropped out trivially. This reflects simply that the scalar product k ■ r is linear in the relative solid angle f2 between 
the vectors whereas the scalar product Tr HK is quadratic in the relative diagonalizing matrix U . The corresponding 
radial Laplacean A^ for N = 1 is identical to the Cartesean A. Therefore, the case N = 1 is t oo tri vial to give any 
further insight. Second, we set N — 2 and find straightforwardly from the differential equation (3.19) 



$f (x,k) 



exp 



(jEi + x 2 )(k 1 + k 2 ) \ (/3+1) ( (xi - x 2 )(k 1 - k 2 ) 



(3.21) 



where the vector Bessel function in d dimensions as defined in Eq. ( |2.3| ). This functions appears in the solution, 

because the differences x\ — x 2 and k\ — k 2 directly correspond to the lengths |r| and |fc|. In higher matrix dimensions 

A, this simple correspondence is lost. However, we will see in great detail that the features of the functions G>£ (x, k), 
in particular whether or not explicit solutions can be constructed, are stronger influenced by (3 than b y N. 

Anothe r po int in this context deserves to be underlined. In the vector case, the differential equation ( 2.10) ) and the 
solution (2.6) were constructed for integer dimensions d. However, both equations are also well defined for any real 
and positive d. Similarly, we observe in the matrix c ase t hat the differential equation (3.19) was derived for the cases 
/3 = 1,2,4. However, neither itself nor its solution ( 3.21 ) for N = 2 are confined to these cases f3 = 1,2,4, they are 
valid for any real and positive j3. Thus, the cases j3 = 1, 2,4 which correspond to a matrix model, i.e. to the defining 
integral ( 3. 14j ) of the matrix Bessel functions, are only special cases of a much more general problem, namely finding 
the solutions of the differential equation (3.19) for every integer N and for arbitrary real values of f3. We will return 
to this in Sec. 0. 
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IV. RECURSION FORMULA 



The matrix Bessel functions show a rec ursive structure which we constr uct by introducing radial Gelfand-Tzetlin 
coordinates. The result is stated in Sec. IV A and derived in Sec. IV B. The corresponding invariant measure is 



calculated in Sec. IV C| 



A. Statement of the Result 



The matrix Bessel functions, defined in Eq. (3.15) 



<S> ( N ] (x,k) = J dfi(U) exp(iTr[/W/c) (4.1) 

depend on the radial space of the eigenvalues x and k. As before, we write x — diag (xx, . . . , Xpf) and k = 
diag (kx, . . . , kpf) for j3 = 1, 2 and, for j3 = 4, we write x = diag (xx, x%, . . . , xn, xn) and k — diag (kx, kx,..., k^, kjsr). 
We emphasize that the radial spa ces do not lie in the manifolds covered by the groups U(N;(3). However, we will 
show that the group integral ( |4.l| ) can be exactly mapped onto a recursive structure which acts exclusively in the 
radial space. This remarkable feature is the main result of this section. 

Under rather general circumstances, the matrix Bessel functions &^\x,k) can be calculated iteratively by the 
explicit recursion formula 

^\x,k) = J dfi(x',x) exp(i(Trx- Trx')k N ) Q^l^x' ,k) (4.2) 

where ^^_ 1 (x' ,k) is the group integral ( |4.1j ) over U(N — l;/3). We have introduced the diagonal matrix k = 
diag (kx, ■ ■ ■ , fcjv-i) f° r P — h 2 and k — diag (kx, kx, ■ ■ ■ , fcjv-i, fcjy-i) for /3 = 4 such that k — diag (k, kjy) for (3 = 1, 2 
and k = diag (k, fc/v, fczv) for [3 = 4. Importantly, the N — 1 integration variables x' nl n = 1, . . . , N — 1, ordered in the 
diagonal matrix = diag (x'-y, . . . , a;^ v _ 1 ) for /3 = 1, 2 and a;' = diag (x' 1; . . . , ar^ v _ 1 , 2:jv_!) for /3 = 4 are arguments 

of $^^(2:', fe). Moreover, we notice that their further appearance in the exponential is a simple one due to the trace. 

The coordiriates x' are constructed in the spirit of, but they are different from, the Gelfand-Tzetlin coordinates 
of RefsJl3'E£l. To clearly distinguish these two sets of coordinates from each other, we refer to the latter as angular 
Gelfand-Tzetlin coordinates and to the variables x' as radial Gelfand-Tzetlin coordinates. The difference is at first 
sight minor, but of crucial importance. In the angular case, x is in the Cartan subalgebra belonging to U(N; [3). In 
the radial case, however, x is in the radial space of the eigenvalues of the real-symmetric, Hermitean or Hermitean 



self-dual matrix H, which are the arguments of the functions (4.1). While the angular Gelfand-Tzetlin coordinates 
never leave the group space, the radial ones establish an exact and unique relation between the group and the radial 
space. The radial Gelfand-Tzetlin coordinates re-parametrize the sphere that is described by the N th column Un of 
the matrix U G U (N; [3). The recursion formula ( fh^ ) can only be constructed in the radial coordinates x', but not in 
the angular ones. The radial and the angular Gelfand-Tzetlin coordinates are, in general, different. They happen to 
coincide for (3 = 2, i.e. for the unitary group U(N). This illustrates, in the framework of our recursion formula, the 
special role played by the unitary group. 

The invariant measure dfj,(x' , x) is, apart from phase angles, the invariant measure <i/i(J7jv) on the sphere in question, 
expressed in the radial coordinates x'. It only contains algebraic functions and reads explicitly 

= (- n<*» - <o) ( " 2)72 . ( 4.3) 

\ x ) \ n,m / 

The normalization constant obtains from results in Gilmore's book0. It ensures normalization to unity according to 



Eq. (3.5). The domain of integration is compact and given by 

x n < x' n < x n+1 , n = 1,...,(N- 1) , (4.4) 

reflecting a "betweenne ss c ondition" for the radial Gelfand-Tzetlin coordinates. This is why no absolute value signs 
appear in the measure fl4.3| ). 
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The general recursion formula (4.2) states an iterative way for constructing the matrix Bessel function <&ffl(x, k) 

for arbitrary N from the matrix Bessel function $2 \ x ,k) for N = 2 which can usually be obtained trivially. We 
remark that the recursion formula allows one to express the matrix Bessel functions in the form 



*%\x,k) = 



N-l 



n=l 

exp 



(n) T (n-1)) 

fif&a^ 1 ) - Tra; (,l) )fcjv_ n+ i) exp (m^^i) 



where we have introduced the radial Gelfand-Tzetlin coordinates it' , 
1, . . . , (N - 1). We define x^ = x and = x' . 



(4.5) 

n on N — 1 levels n = 



B. Derivation 



We introduce a matrix V = diag {V, V ) with V" G U(N-1; j3) and y G 17(1; j9) such that V £ U(N-1; P)®U(1; /?) C 
[/(TV; /3) and multiply the right hand side of the definition (4.1) with 



1 = / dfi(v) = j dfi(v ) J duty) . 

The invariance of the Haar measure dfi(U) allows us to replace U with UV' and to write 
<P { x\x,k) = j dn{V ) j d(j,(V) j dfi(U) expiiTrU^xUV^kV) . 



(4.6) 



(4.7) 



We collect the first N — 1 columns U n of U in the N x (N — 1) rectangular matrix B such that B = [Ui U2 • • ■ Utt—i] 
and U = [B Un]- We notice that 



B*B = 1 

A 

BB^ = U n Ul = 1 N - U N U N 



LJV-1 
N-l 



(4.8) 



As already stated in Sec. |III A| , the elements of a vector or a matrix are scalar for (3 = 1, 2 and quaternion for (3 = 4. In 
this sense, we also write In as the unit matrix for f3 — 4 because its elements are r 1 - - 1 . By defining the (N — l) x (iV— 1) 



square matrices H = B^xB and if = V^kV we may rewrite the trace in Eq. (4.7) as 

Tr U^xUV^kV = Tr HK + H NN k N 



(4.9) 



with Hnn = U\j xUm according to Eq. (3.6). We notice that Vq has dropped out. Since the first term of the right 
hand side of Eq. (4.9) depends only on the first N — l columns U n collected in B and the second term depends only 
on Un, we use the decomposition 



dfi{U) = dn{B)d(i{U N ) 
of the measure to cast Eq. Q4.7D into the form 

/ dfi(Uff) exp{iH N Nk N ) / dfi(V) / dfi(B) exp(«Tr/j£) 



(4.10) 



(4.11) 



where we have already done the trivial integration over Vq. 

The difficulty to overcome lies in the decomposition (4.1C). While dfi(UN) is simply the invariant measure on the 
sphere described by Un, the measure dfi(B) is rather complicated. Pictorially speaking, the degrees of freedom in 
dfi(B) have always to know that they are locally orthogonal to Un- Thus, dfi(B) depends on Un- Luckily, there is 
one distinct set of coordinates that is perfectly suited to this situation. It is the system of the radial Gelfand-Tzetlin 
coordinates. We construct it by transferring the methods of Ref.o for the angular case to the radial case. 
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The N x N matrix (ljv — UnII n ) is a projector onto the [N — 1) x (AT — 1) space obtained from the original N x N 
space by slicing off the vector Un. We project the radial coordinates x onto this space and study its spectrum. The 
defining equation reads 

(ljv - U N U N ) x (1 N - U N U N ) E' n = x' n E' n n = 1, . . . , N - 1 . (4.12) 

Equation ( f,12| ) determines the N — 1 radial Gelfand-Tzetlin coordinates x' n and the corresponding vectors E' n as 
eigenvalues and eigenvectors of the matrix (ljy — UnU n ) x (In — UnU n ) which has the rank N — 1. Since we have 
by construction U N E' n = 0, we may as well write 

[l N -U N U N )xE' n = x' n E' n , n=l,...,7V-l . (4.13) 
The eigenvalues x' nl n = 1, . . . , N are obtained from the characteristic equation 

= Dot ((ljv - U N U N )x - 
= Det (x — x' n ) det ^1^ — (x — x' ji )~ 1 UnU n x^ 

= Det (x - x' n ) ( 1 - U N ^ T U N ) 

= -x' n Det (z - a^) Tr U^-^Un . (4.14) 

a; x n 

Together with the normalization Tr U n Un = 1, this yields the N equations 

N 0-1 



l = TrU N U N = 



(a)2 
n— 1 a=0 



AT /3-1 (a)2 

= T4^ = EE^7. n = l ) ...,AT-l. (4.15) 



m— 1 a— 



In these formulae, the trace Tr is only needed in the symplectic case. We notice that the equations for the variables 
x' depend on the variables x as parameters. We emphasize once more that x in these equations is in the radial space 
and, in general, not in the Cartan subalgebra of Z7(iV; (3). 

At this point, it is not clear yet why the introduction of the radial Gelfand-Tzetlin coordinates is at all helpful. 
The great advantage will reveal itself when we express the matrix H and the matrix element Hnn in the trace (4.11) 
in these coordinates. To this end, we first multiply Eq. (4.12) from the right with E') t and sum over n, 

N-l 

(l N -U N U N )x(l N -U N U N ) = £y n £££# (4.16) 

n=l 

where we used the completeness relation 

N-l 

]T E' n E% + U N U N = 1 N . (4.17) 



Taking the trace of the spectral expansion (4.16) we find immediately 

Tra: - Trx' = TrU N xU N = H NN (4.18) 

This is a remarkably simple result. An analogous expression exists for the. N . N matrix element of the unitary group 
in the theory of angular Gelfand-Tzetlin coordinates for the unitary groupllaEa. Here we have shown that Eq. ( |4.1§| ) 
is a general feature in every radial space. 

We now turn to the (N — 1) x (N — 1) matrix H . Its N — 1 eigenvalues y n , n = 1, . . . , N — 1 arc determined by 
the characteristic equation 
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Det (H - y n ) = Det (B^xB - y n 



Vn 
1 

Vn 



Det (BB^x - y v 



Det 



nliv - UnU^x - y„) 



(4.19) 



where we us ed Eq. (4.8) and re-expressed a (N — 1) x (N — 1) determinant as a N x N determinant. The comparison 
of Eq. ( 4.19| ) with Eq. (4.14) shows that, most advantageously, we have y n = x' n , n — 1, . . . , N — 1. Thus we may 
write 



H 



B^xB 



U^x'U 



(4.20) 



by introducing the (N — 1) x (N — 1) square matrix U which diagonalizes H. Obviously, U must be a complicated 
function of the N x (N — 1) rectangular matrix B, i.e. of the columns U n , n = 1, . . . , N — 1. However, all we need to 
know is that U must be in the group U (N — l;/3) because, by construction, H has the symmetry = H. 
Collecting everything, we cast Eq. (4.11) into the form 



^\x,k) = J dfi(x',x) exp{i(Trx-Trx r )k N ) 

dn{V) / dfi(B) exp(iTrU^x'UV 1 kV) 



(4.21) 



We may now use the invariance of the Haar measure d/i(V) to absorb U such that 

^\x,k)= I d(J,{x',x) exp(i(Trx-Trx')k N ) 

dfi(V) exp^Trx'V^fcVO / d^B) . (4.22) 

Thus, the integration over B is trivial and yields unity due to our normalization. The remaining integration over 
V gives precisely the matrix Bessel function ^^^(a;', fc). This completes the derivation of the recursion formula 
in Sec. [IV A . The reader experienced with group integration has realized that the introduction of the matrix V = 
diag (V, Vq) was not strictly necessary. Alternatively, one could have shown that the measure d/i(B) can be identified 
with dfi(U) and have done the corresponding integral. However, we believe that the introduction of V makes this 
part of the derivation more transparent. 



C. Invariant Measure 



The invariant measure d^U^) has to be expressed in terms of the radial coordinates x' . To this end, we first have 
to solve Eq. (4.15) for the moduli squared of the vector Un as a function of the new coordinates x' . Since Eq. ( 4.15| ) 
for the total moduli square for all (3 coincides with-the equation for the angular Gelfand-Tzetlin coordinates of the 
unitary group, we can use the results as derived inErE3. We have in the three cases 



\U n 



0-1 r\ N ~ 1 ( - 1 

a=0 



Y[m^n( Xn Xm ) 



(4.23) 



The betweenness condition (4.4) follows from the positive definiteness of this expression. We parametrize the remaining 



degrees of freedom of U u n in the cases j3 = 2,4. We set U^ N = cos7„ and U^ N = sin7„ in the case 13 — 2 and 



U, 



N 



cos tjj n exp(i7i 1) ) sin ip n exp(i7^ J ) 



(2h 



sin^n exp(-i7i 2) ) cos ip n eTq>(-ijn>) 



(4.24) 



for (3 = 4 in the basis (3.2). The invariant length element reads 
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N 0-1 



Tr dUjfdUjf = EE 



n—l a—0 

N 



r(«) 
J nN 



0/2 



n—l ^ ' ' z—1 



(0^2 



+<5 /34 |[/ nJV | 2 (dcos</> I1 ) 2 J . (4.25) 
To express the differential d\U n N\ 2 in terms of the dx' n , we again take advantage of the results in Refs.0El 

N , N-l T-rN-l, I _ I \ 

^ 4\U nN \* ^ 4ll„ = i(^ " <) 

From these equations, we can read off the metric g in the basis of the coordinates x' n , 7„ and ip n . Conveniently, it is 
diagonal. The determinant of g is given by 



det g — 



(4.27) 



which yields the invariant measure d/i([/jv) in terms o f the x' n and of the additional coordinates 7^ and ip n . These 
angles can be integrated out trivially. This yields Eq. (4.3). 



V. RADIAL FUNCTIONS FOR ARBITRARY (i 



Remarkably, the recursion introduced in the previous section, is the iterative solution of the radial equation for 
arbitrary values of j3. Thus, the matrix Bessel functions are special cases of more gener al functions whic h we want to 
refer to as radial functions. We give the precise fo rmul ation of the problem in Sec. V A and show in Sec. VB that the 
recursion is the general iterative solution. In Sec. V C , we discuss a Hankel ansatz for the radial functions. 



A. Definition by the Differential Equation 



In Sec. IIIB , we defined the matrix Bessel function through the group integral (3.14) or, equivalently, the group 
integral ( 3.1 5| ) . This definition confines the dimension f3 to the values (3 = 1, 2, 4, corresponding to the groups U (N; (3). 



However, discussing the simplest case N = 2, we already saw in Sec. |III B| that $2 is wen defined for arbitrary values 



of (3, This was a simple consequence of the explicit form (|3.21 ) which expresses $ 2 m terms of the Bessel function 



X 



(0+1) 



The latter is known to be well defined for arbitrary (3. Hence, we conclude that the cases [3 — 1,2,4 which 
relate to matrices and groups are embedded into a space of far more general functions. 

It seems natural that this phenomenon also extends to N > 2. The problem has to be posed as follows: We seek 
the solutions <I>^ of the differential equation 



N 



-E fc 



2 k) 



(5.1) 



where the operator is given by 



A x = 



N B 2 

T — 



' i"- — r r— ' dx n 



dx r , 



(5.2) 



Here, f3 is arbitrary. For technical reasons, however, we restrict ourselves for the time being to real and positive values 
of (3. We make no reference whatsoever to matrices, eigenvalues and groups. To emphasize this, we view x and k as 
sets of N variables x n , n = 1, . . . ,N and fc„, n = 1, . . . ,N for every positive (3. We do not use traces. 
We require that the solutions are symmetric in the argument 
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*%>(x,k) = < J (fc,x) (5.3) 

and normalized 

$^(0,fc) = 1 and *^(ar,0) = 1 (5.4) 

at the origin x — and k = 0. 

In the sequel, we want to refer to the functions <& £ (x, k) for arbitrary (3 as radial functions while we reserve the 
term matrix Bessel functions to the cases f3 = 1,2,4 where the direct connection to matrices and Lie groups exists. 



B. Recursive Solution 



We claim that the solutions are, for arbitrary (3, given as an iteration in N by the recursion formula 

^\x,k) = / dfx(x',x) exp ( i I >> - >>' | few | *%Utf,k) (5.5) 



J dfj,(x',x) exp ^i (y^^ x ~ X] x '^j kl ^j ®N-i( x '' k ) 



where ^>^_ 1 (x' , k) is the solution of the differential equation (5.2) for N — 1. Here, k denotes the set of variables 
k n , n = 1, . . . , (JV — 1) and x' the set of integration variables x^, n — 1, . . . , (N — 1). The integration measure 



A^(x) 




is the continuation of Eq. (4.3) to arbitrary positive (3. The normalization constant 



G W = 2 n-i W/2) 



r JV (/?/2) 



(5.6) 



(5.7) 



is also the continuation of the constant in Eq. (p~3|). We calculate it in App. |]. As in the cases (3 — 1,2,4, the 
inequalities 



n+l 



n=l,...,(A-l) 



(5.8) 



define the domain of integration. 

We stress that we derived the recursion formula (5.5) in Sec. IV for the cases (3 = 1,2, 4. To prove that it is the 
iterative solution for arbitrary positive (3, we show that it solves the differential equation (5.1). The keystone for the 
proof is the identity 



A x ^\x,k) 



-k%^\x,k) 



+ J d/i(x',x) exp I i y 

A^l^x^k) 



N N-l \ 

Xn ~ X/ X ™ ) 

n— 1 n— 1 / 



(5.9) 



which is derived in App. Equation (|5.9| ) establishes a not immediately obvious, but nevertheless natural connection 
between, on the one hand, the actio n of the Laplacean in the N variables x n on the radial function in N dimensions, 
i.e. on the recursion integral (5.5), and, on the other hand, the recursion integral over the Laplacean A x i in the 
N — 1 variables x' n acting on the radial function in N — 1 dimensions. There is a compensation term which is just 



kl^MxJt). Thus, we can prove the eigenvalue equation (pTJ) by induction: assuming that it is corr ect for N — 1, 



implies Eq. (5.1) for N. The induction starts with N — 2 wher e the eigenvalue equation (5.1) is clearly 



identity (I 

valid for arbitrary (3 as sho wn in Sec. MB by deriving the explicit solution ( 3.21 

The sym metry relat ion (5.3) is non-trivial. In the matrix cases [3 = 1,2,4, it is obvious from th e in tegral defini- 
tions (3.14) and (3.15). For arbitrary f3, we cannot use this argument, we only have the recursion (5.5). In App. ^, 
we prove the symmetry relation (|5.3|) by an explicit change of variables. 



The normalization $^)(x, 0) = 1 in Eq. (5.4) follows directly from the normalization of the measure (5.6). The 



symmetry relation (5.3) then also yields <f>^ (0, k) = 1. 
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Regarding the domain of /3, a comment is in order. We have seen in Sec. Ill B| that for N — 2 th e m atrix Bessel 
function is well defined for arbitrary complex (3. This should also be true for our recur sion formula (5.5). However, 
for (3 < non-integrable singularities arise at the boundaries in the integral in Eq. (5.5). At the same time the 
normalization constant becomes zero for (3 = 0, —2, —4, . . . compensating the singularities of the integral. This makes 
the recursion formula for (3 < not ill-defined but it gets more difficult to treat. Therefore, we have restricted 
ourselves to positive values of (3. . 

In the work of Okounkov and 01shanskiE3 an expansion of the radial functions for arbitrary (3 in Jack polynomials 
is derived. The series run over sets of partitions {(J.}. These authors also derive a recursion formula for the Jack 
polynomials depending on one set of continous variables x, say, and belonging to such partitions {/i}. It is related 
to, but different from ours which involves two sets of continous varables x and k. The crucial difference rests in the 
exponential function which is present in our formula (5.5), but n ot in the formula of Rcf.EJ. Importantly, it is this 
exponential term which makes sure that the symmetry condition (5.3) is fulfilled on all levels of the recursion. Since 
the Jack polynomials themselves do not obey such a symmetry condition, there is no exponential term in the recursion 
formula of Ref.e9. However, it must be possible to derive the recursion formula for the radial functions from the one 
for the Jack polynomials. An interesting, although probably not very elegant papproach would be the following: If 
one inserted the recursion formula for the Jack polynomials into the expansionE] of the radial functions in terms of 
these Jack polynomials, one ought to see that the series over the partitions can, at least partly, be resummed to yield 
the exponential function present in the recursion formula (5.!:). This is remarkable and could be very helpful for 
the application of Jack polynomials, because, in general, resummations over partitions are known to be difficult and 
involved. For the connection to Calogero-Suthcrland models, we refer the reader to Sec. VII. 



C. Hankel Ansatz 



In the spirit of Eq. (2.12) for the vector case, we make a Hankel ansatz for our radial functions for arbitrary positive 
(3. We also do this i n vi ew of the applications in Sec. VI. Since the sum over the on the right hand side of the 
eigenvalue equation (pTJ) is invariant under all permutations of the k n or, equivalently, their indices n, we can label a 



set of solutions ^^' u (x, k) by an element u) of the permutation group Sn of N objects. For these solutions, we make 
the ansatz 



exp (iJ2n=l x nku,( n ) 

\A N (x)A N (k)\P/z 



Wghx,k). 



(5.10) 



where uj(k) is the diagonal matrix constructed from k by permuting the k n , or the indices n. The full solution 



(x, k), satisfying the constraints (5.3) and (5.4), is then, apart from possible normalization constants, given as 
the linear combination 



(5.11) 



of the functions ( 5. 10] ). Here, tt(oj) is the parity of the permutation. 



We find for the function W^\{x, k) the differential equation 



L wMk) W<ghx,k) = 



(5.12) 



where the operator is given by 



L, 



N & 
T — 

71=1 " 



n— 1 x / n<rn K J 



(5.13) 



This differential equation generalizes Eq. ( 2.13 ) to the matrix case for (3 = 1,2,4 and, furthermore, the latter to 
general radial function for arbitr ary (3. 



Again, due to the symmetry (5^), the differential equation ( 5.12 ) must also h old if x and u>(k) are interchanged. 
It is the last term of the operator L x u ^ that makes the differential equation (5.12) so difficult. This shows that 
the case (3 = 2 corresponding to unitary matrices U S U (N) is special: the last term vanishes and we simply have 
^ / ivL( x ' = 1- This is the Itzykson-Zuber caseil. For arbitrary (3, it is obvious from the differential operator that 
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Wjfl(x,k) -» 1 if 1^ 



oo for all pairs n < m. Once more, this must also be true if \k n 



Thus, 



we expect that Wj^\,{x, k) is some kind of asymptotic series, generalizing Eq. ( |2.14|) in the vector case. 

Hence, we re- derive a known result by concluding that the leading contribution in an asymptotic expansion of the 
functions (5.10) is given by 



exp 



(n) 



\A N (x)A N (k)\P/i 



According to Eq. (5.11), this means that 



^\x,k) 



det [exp(ix n k 



m ! \n.m- 



\A N (x)A N (k)\W 



and I 



(5.14) 



(5.15) 



are large for all 



is the asymptotic behavior of the radial functions (x, k) if the differences \x n — x Tt 
pairs n < m. 

The functions W$l(x, k) are translation invariant, i.e. they depend only on the differences (x n —x m ). We show this 

in App. Due to the symmetry, this argument carries also over to k and W^\{x, k) depends only on the differences 
(fc n — k m ) as well. Moreover, the symmetry implies that it depends only on the products (fc w („) — k u ^ m -\)(x n — x m ). 
Collecting all these pieces of information, we make the ansatz 



w { N p hx,k) 



E 



a 



A»12A 1 13'"M(JV-1)1V 



{ M } Iln<m ((^w(n) ^u{m))( x n x m)) 



(5.16) 



with coefficients o J ^ 12tll3 ... tl{N _ 1)N that depend on N(N — l)/2 integer indices fJ, nm , as many as there are differences. 
The summa tion is over the set of these indices. The presence of the k n makes it very difficu lt to solve Eq. (5.12) with 
the ansatz ( 5.16 ). In the vector case, one easily sees that the differential equation ( 2.1 3| ) in r can be transformed 
into an eq uatio n in the dimensionless variables kr such that k does not appear anymore. This leads to the simple 
recursion (2.15) for the coefficients. Here, in the matrix case, the k n cannot easily be absorbed and the recursion 
formulae for the coefficients will depend on the k n in a non-trivial way. However, in some simple cases, it is possible 
to solve them. These difficulties were an important motivation for us to develop the methods which we introduced in 
Sec. Evl 



VI. APPLICATIONS 



Can we obtain explicit formula for the radial functions by using the recursion formula (5.5) ? — At least in some 
cases, this ought to be possible. Here, we present our first attempts. 

For the sake of completeness, we comment once more on the special case (3 — 2, i.e. the unitary case. Obviously, 



the measure (5.6) simplifies enormously. This is so because the radial Gelfand-Tzetlin coordinates coincide with.ihe 
angular ones. Thus, the case (3 = 2 is identical to the re-derivation of the Itzykson-Zuber integral by ShatashvilO. 
We now consider the orthogonal case 0=1. The recursion formula reads 

J d~ lln,m( Xn ~ X 'm) 

/ / N N-l \ \ 

ar l k N j *$_ 1 (* , ,fc)dM • (6.1) 

The square roots appearing in the measure make a further evaluation very difficult. As obvious from the trivial case 
N = 2, given in Eq. (3.21), the function $^(a;, k) will be an infinite series for all values of N. However, we expect 




that, due to the different construction, this series is different from the expansion in zonal functions which was obtained 
by MuirheacM 

Obviously, there is a pattern emerging. The integration measure ( [5.6| ) is purely rational for all even and positive 
values of (3. This is reminiscent of the situation for vector Bessel functions in odd dimensions d, which consist of a 
finite number of terms, as discussed in Sec. [□]. Hence, we conjecture that the radial functions $^' ) (x, k) can also be 
written as a finite sum, exclusively containing exponential and rational functions. 
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For all other values of/3, the measure ( |5 ,6| ) is algebraic, but not rational, and the radial functions <&jy (x, k) must be 
infinite series. Nevertheless, these infinite series contain exponential and rational functions. Thus, they are different 
from expansions in terms of zonal polynomials. 

To furnish our conjecture about the form of the radial functions for even and positive values of (3 with an illustrative 
example, we turn to the unitary-symplectic case /3 = 4. To simplify the notation we avoid the imaginary unit by 
writing 

$$(-ix,k) = / exp(Tru -1 a:uifc)<^(Z7) , (6.2) 

JUeUSp(2N) 

where x = diag (xi, X\, . . . , xn, xn) and k = diag (ki, k\, . . . , fcjv, k^) are diagonal matrices with Kramers degeneracies. 
The starting point of the recursion is the smallest non-trivial case N = 2, i.e. the group USp(4). We obtain after an 
elementary calculation 



^\-ix,k) 



- 2 

G 2 4) E ( A 2 M A2 ~ A3 M A3 I^mS ) exp(Trxc(fc)) . (6.3) 



Ai(x)AlKfc)) A3 (x) ( w (fc)) 



The sum runs over the elements of the permutation group Sat for N — 2. Inserting Eq. (6.3) into the recursion 
formula, we find for USp(6), the next step in the recursion, 



2 i r 

Xa 



exp 



Al(x)Al(Lo(k)) 
2 



f (Trx — Trx')k 3 + Tr xuj(k)] f - — — — - — 

V 'V I A 2 (x>) A|(x')A 2 (u;(fe)) 



(6.4) 



Although the integrand is finite everywhere, in particular at the denominators A2(x') and A|(x') raise 
a technical difficulty. The key to remove them is to use the identity 

-*-=-(-*- -°- )^ (65) 

A|(x') \dx[ Ox'J A 2 (x') { ' 

and to observe that the product Yii=i Tij=i( x i ~ x j) annihilates all boundary terms. Hence, we can integrate by parts 
and arrive at 

^\-ix,k)=G^j: 



weS2 Al(x)AUu>(k)) 

dxi / gLx 2 ^ (xj - x\ ) (xj - x' 2 ) 

-1 i=1 3 = i 

((Trx-Trx')fc 3 + Trxw(fc)) , (6.6) 



where no denominator is left. Due to the permutation symmetry of the original int egral, we can restrict ourselves to 
the unity element e of the permutation group in the further evaluation of Eq. (6.6). Thus we need only to consider 
the limits x\ — > X{, i = 1, 2 while integrating by parts. After collecting orders in k we find 

d> 3i£ ( lX,k) - G 3 G 2 A 3 (x)A 3 (fc) 




3 

(xi — Xj)(ki — kj) + 12 | exp (Trxfc) . (6-7) 

i<j 
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By introducing the composite variables 

= {xi - Xj)(k u (i) - k u (j)) i,i = 1, . . - , 3, uj € S 3 , (6.8) 



we can express $3 (— ix, k) in a compact form as 



*W(-fa,fe) = G«<# ) £ 



1 



^ A|(x)Al( W (*)) 
3 



4 + [J(2 - cxp (Trxw(fc)) . (6.9) 



So far, we have not been able to extend this procedure to all values of N. 

However, we succeeded in calculating ix, k), i.e. the case of the group USp(8), by an hybrid method which 

combines informations obtained from the recursion with an Hankel ansatz as described in the previous section. We 
extend the right hand side of Eq. (|6.9|) for N = 3 to N = 4 and use this expression as an ansatz for the function 



W^ u (x, k). As it turns out, a correction term is needed and, furthermore a correction to the correction. Fortunately, 
there is a structure to this. We give the details in App. ||. We emphasize that the knowledge of $ 3 4 ' ) (— ix, k) 
is essential for this hybrid procedure, in particular the fact, that $g 4 ' (— ix, k) contains only linear terms in every 
composite variable z u ^j\. Up to a normalization, (— ix, k) is given by 

*i 4) H^) = E A3 W AKo;( fc )) (n( 2 -Mii)) + 

\ e n (2-m«))+ 

Km<n i<3 

lE II (2-^ fo) ))exp(Tr a ;c i ;(fc)) . (6.10) 

Km i<j,^lk,^ln J 

Comparing this result with Eq. ([T{]) we notice that, once more, the composite variables enter only linearly 

in the polynomial part of <frf\~ix, k). Similarly, the spherical Bessel function ji(z), which is the counterpart of 



<!>$(x, k) in the vector case given in Eqs. (2.6) and (3.21), has a polynomial part linear in z. We expect that such 



analogies are also present for higher values of (3 and the dimension d. 

Formula ( 6.10| ) indicates a general structure for <&j^(x, k). The leading term is always the generating function of 



the elementary symmetric functions in z. To this term combinations of other symmetric functions are added, where 
certain combinations of indices are cut out. 

VII. CONNECTION TO CALOGERO-SUTHERLAND MODELS 

The radial functions are related to, but different from, the eigenfunctions which are usually employed in models of 
the Calogero-Sutherland type. Since this issue is so important for applications and so often raised in discussions, we 
briefly collect the main points. 

The radial Laplace operator defined in Eq. (|5.2| ) is closely related to the Calogero-Suthepland Hamiltonians. In 
general one can always cast a Fokker-Planck operator in a Hamilton operator by adjunctiona with the square root 
of the stationary probability distribution defined through A x P eq (x) — 0. Choosing P e q(^) = |Aat(:z;)|~ /3 / 2 , the 



operator (5.2) can be associated with the Hamiltonian 



N ^2 n,n ^ N 



02 p { p_ 2) » 1 

It describes a scattering system with a continuous spectrum, the large time behavior is determined by the states near 



the ground state. Apart from a sign, this operator Hd coincides with the operator L x ,o in Eq. (5.13) for k = 0. We 
also notice that the interaction vanishes for (3 = 2. 
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To have a well defined thermodynamic limit one often confines the motion of the particles to a circle. This yields 
the Calogero Sutherland Hamiltonian 



N q2 «n N 



H V ° 2 I M ~ 2) V {7r/N)2 (72) 

which can also be derived directly from Dyson's circular ensembles0. Another way of confining the particles is by a 
harmonic potential. This leads to the Calogero HamiltoniarJa 



-1 ft n<m v ' 



In the thermodynamic limit the particle density Ri(x) of the ground state is described by Wigner's semi-circle law. 
The mean particle level spacing D — l/i?i(0) scales as D oc l/V^V. Therefore in the thermodynamic limit the 



harmonic confining term in Eq. (7.3) vanishes on the scale of the mean level spacing. On this unfolded scale the 
correlation functions become independent of the confinement mechanism. The three Hamiltonians He, Hp and Hcs 
are known to be integrable systems for arbitrary /£3. However, the three values f3 = 1,2,4 are distinguished, since 
they establish a connection to the random matrix ensembles. Indeed, for these values of (3 they belong to a much wider 
class of integrable systems, which canJae constructed by means of the root space of a simple Lie algebra or - still more 
generally - of a Kac-Moody algebraic. This class comprises Hamiltonians which can be derived by an adjunction 
procedure from a Laplace-Beltrami operator of a gWHp acting in a symmetric space. This space has positive curvature 
for Hqs and zero curvature for Hp, He- In Refs.oa it was pointed out, that the Dorokhov-Mello-Pereyra-Kumar 
equation for scattering matrices with broken time reversal symmetry corresponds to a Laplace-Beltrami operator in 
a symmetric space of negative curvature. 

Eigenfunctions iPn e(x) of the Hamiltonians He, Hp, Hcs with eigenenergy E for arbitrary (3 are known. Essen- 
tially these solutions are products of the ground state wave function and symmetric polynomials in the coordinates 
x of the-|iY| particles. In case of the Calogero-Sutherland Hamiltonian Hcs, these polynomials are the Jack polyno- 
mialsE§BO. In this approach, the energy eigenvalues E are labeled by a partition of length N. The crucial difference 
to the matrix Bessel functions is that the Jack polynomials are symmetric polynomials in one set of variables x only 
whereas the matrix Bessel functions $^(x, k) are symmetric in two sets of variables x and k. Importantly, they 

are, in addition, symmetric under interchange of the two sets of variables, <&^'(x, k) — $^\k,x). This is reflected 
in the fact that the operator L x>UJ (k) emerging in the Hankel ansatz depends on k while Hjj does not. Due to their 

symmetry, the matrix Bessel functions <5>^\x, k) are, at least for Hjj, the more natural eigenfunctions. This becomes 

obvious in the fact that, to obtain orthogonality conditions, one has to sum the i/^e^) over an infinite number of 

partitions. On the other hand, orthogonality relations are an inherent feature of the <fr^\x, k) due to their meaning 
in the Fourier-Bessel analysis, as discussed in App. |C| 

In other words, the functions can be viewed as a basis in an expansion of the 3>^\x, k). We can consider 

the variables k as a set of real numbers corresponding to the energies E labeling the eigenstates of Hp. The matrix 



Bessel functions (3.15) arc solutions of the Schrodinger equation with Hamiltonian Ho for the coupling parameters 



[3 = 1,2,4. The recursion formula (5.5) represents an analytic continuation of these integral solutions to arbitrary 

,(/3) 
N 

no infinite restimmation is required to obtain functions of the type ^^\x,k). Nevertheless Forresteri and Nagao 



positive p. All these functions (x, k) have, for arbitrary (3 additional features, such as the symmetry in x and k, 
which have no analogue in the functions tp^ E (x). The merit of our recursion formula lies in the fact that, a priori, 



and Forrester^ showed that such resaimmed expressions can successfully be used in certain cases. They treated the 
case of Poissonian initial conditionsEE for the Calogero-Sutherland Hamiltonian Hcs and derived exact expressions, 
for the correlation functions for arbitrary (3 for one or two particles. This is also related to the works of MuirheadEj 
and PandeyEil 



VIII. SUMMARY AND CONCLUSION 



We presented a recursive construction for certain spherical functions. We referred to them as matrix Bessel functions 
because, first, they are a natural extension of vector Bessel functions in the sense that the integration over a group 
corresponds to the integration over a solid angle and, second, they satisfy a partial differential equation generalizing 
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the Bessel ordinary differential equation. For matrices, the index j3 labeling the groups appears analogously to the 
dimension d in the case of vectors. The introduction of radial Gelfand-Tzetlin coordinates, which are related to but 
different from the ordinary angular ones, was crucial for the recursion. The Cayley transformation ought to provide a 
connection between the angular and the radial Gelfand-Tzetlin coordinates^ As evident from its construction, the 
recursion maps an integral over a group fully onto an iteration which exclusively takes place in the radial space. 

Remarkably, the recursion turned out to be far more general than was to be expected, at first sight, from the proof 
which involved Lie groups. We showed that our recursion is also the iterative solution of the corresponding partial 
differential equation for arbitrary values of (3. We introduced the term radial functions for this generalization of 
matrix Bessel functions. We expect that one has to employ the theory of quantum groups to give a group theoretical 
derivation of the recursion formula for arbitrary values of f3. 

Using the recursion formula, we discussed the structure of radial Bessel functions. We conjectured that, for even (3, 
they can be written as finite sums involving only exponential and rational functions. We illustrated that by working 
out, for (3 — 4, the cases of N — 3 and N = 4 distinct eigenvalues. Further evaluation of explicit formulae for arbitrary 
N and, maybe, for all even (3 does not seem impossible. Work is in progress. The extension of the stationary phase 
approach by Duistermaat and HeckmanO to higher orders could, for even (3, be an alternative to derive such explicit 
results, because the expansion terminates. In this context, we mention that the radial function for higher values of 
(3 are, to some extent, but not fully, the higher order radial functions for lower values of (3. This also generalizes the 
situation for ordinary Bessel functions. However, there are many more higher order radial functions, they are not at 
all exhausted by this mapping between values of (3. 

In the present contribution, we only focussed on ordinary spaces, i.e. spaces which are built upon commuting 
numbers. In a second studylij we also address superspaces which involve commuting and anticommuting variables. 
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APPENDIX A: ALTERNATIVE INTEGRAL REPRESENTATION FOR VECTOR BESSEL FUNCTIONS 



The Bessel functions x^^(fcr) arc defined as an integral over angles in Eq. (2.1), but they can also be written as 
integrals over the entire or half real axisa. To make possible an instructive comparison with the matrix case, we quote 
and re-derive the representation 

(d)n , r(d/2) / 2 \ ( ^ 2)/2 / kr f 1\\ dt , A „, 

The singularities have to be treated properly. 

The position vector in the d dimensional space is f — re r where the vector e r parametrizes the unit sphere. To 
integrate over its orientation, i.e. over the solid angle Q, one can re-express the measure as 

dn = ^ff- 5 {el - 1) d d e r . (A.2) 

Here, the vector e r is re-interpreted: its components live on the entire real axis and the domain of integration is the 
full d dimensional space with the Cartesean measure d d e r . The 8 distribution confines the vector to the unit sphere. 



Writing this distribution as a Fourier transform, we obtain from Eq. (2.1) 

/W ^_ r ( d /2) 1 



(kr) = J /2 — I dt d d e r exp (it (e£ - l)) exp (ik ■ re r ) 



V(d,2) r + -^Ht)f_ i (krl ]dt 



./_«, (ity/ 2 1 \ At 

where the integral over e r gave a Gaussian in d dimensions. The contour for the integration over t has to be chosen 



appropriately. Upon a trivial change of variables, this result yields Eq. (A.l) 
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APPENDIX B: ALTERNATIVE INTEGRAL REPRESENTATION FOR MATRIX BESSEL FUNCTIONS 



The matrix Bessel functions 3>^\x,k) for (3 — 1,2,4 can be be written in an alternative way. Although w e can 



hardly believe that this representation is completely new, we could not find it in the literature. Similarly to Eq. (A.l) 
in the vector case, we can write 



03), 
N 



(x, k) = [ d[T] exp (iTr T) Det - f3/2 fx ® jfe - T ® ljv) 



(B.l) 



where 1 at is the N x N unit matrix. The normalization constant is given by 



A 



(0) 

N 



i 0N 2 /2 7r /3N(N-l)/4 
pN+/3N(N-l)/2 



N 



n r w 2 ) • 



(B.2) 



71=1 



The matrix T in Eq. (B.l) is real symmetric, Hermitean or Hermitean self-dual, respectively, for f3 — 1,2,4. The 
measure d[T] is Cartesean and given by Eq. ( |3.3| ). All independent variables in T are integrated over the entire real 
axis. To ensure convergence, the diagonal elements of T have to be given a proper imaginary increment. We notice 
that x and T are N x N matrices for j3 = 1,2 and 2N x 2N for (3 = 4 with doubly degenerated eigenvalues. The 
matrix k is, in all three cases 0, jus t the N x N matrix k — diag [ki, k2, ■ ■ ■ , fcjv), as defined following Eq. (3.10). 
The integral representation (BA) leads to an interesting integral equation for the matrix Bessel functions, 



$%\x,k) = B^Det^^x J d[t]\A N (t)f 



IlTi.m {km in) 



0/2 ' 



(B.3) 



where the normalization constant reads 



" x — 



(B.4) 



The t n in Eq. (B.3) have a proper imaginary increment and their domain of integration is the real axis. Due to the 
symmetry relation ( 3.16| ), the variables x and k can be interchanged in Eqs. (B.l) and (B.3). 

It is not difficult to see from the integral equation (B.2) that the product in the denominator of its right hand side 
can, in the case (3 = 2 be written as 



B 



(2) 
N 



det[S(x n - t m )] 



m I m,m= 



.N 



Y\n,m {^m t n ) 



\A N (k)A N (t)\^ 



(B.5) 



For j3 2, the term Det 1 ^ l3 ^ 2 x contributes. Nevertheless, the product still shares fe atur es with a S distribution. 

To derive this alternative integral representation, we proceed analogously to Eq. ( A.2| ) by re-writing the invariant 
measure of U € U (N; (3) using 8 distributions. The invariance simply means that all columns U n , n — 1, . . . , N are 
orthonormal, Tr U}JJ m = 5 nm . The trace Tr is only needed for j3 = 4, because the entries of U are quaternions in this 
case. Thus, we may write 



N 



dn(U) = Mf d[U] J] 5 (Tr U\ U n - l) J[ S (Tr UlU m ) 



(B.6) 



where d\U] is the Cartesean measure of all entries of U and the integration is for all variables over the entire real 

axis. The constant M#> will be determined later. Ulla« used such forms for the measure to work out certain 
probability density functions. The bilinear forms in the 5 distributions have (3 components for n ^ m, 



0-1 



UlU m = H U m] 



(«) 



.(a) 



(B.7) 



a=0 



We notice that [UlU n 



(a) 



for a > in th e ca se n = m, because the length of every vector is real. Thus, because of 



Eq. ( |B.7| ), the 5 distributions in the measure ( B.6 ) have to be products of S distributions for every non-zero component 

[ulu„ 



We now introduce Fourier representations 
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8 ( [UXU m ] (Q) ) = \ J J dT<& exp (-<2 [l£ f/ ro ] (o) T« 
5([C7K](°)_i) =iy' + JdTWexp(-i([yK]( )_i)T( 2) (B.8) 

for n ^ m and n = m, respectively. The Fourier variables form the elements 

0-1 

T nm = J2 T ™ r(a) ■ ( B - 9 ) 

of a matrix T which is real-symmetric, Hermitean or Hermitean self-dual according to f3 = 1,2,4. We notice that 
the diagonal elements T nn = T„„ are always real. 

S(TiUlU m ) =-£ /" e^T^ exp (-iTr Ul(T nm ® l N )U m 

-iTrUl(T: m ®l N )U n ) 

1 /■ +0 ° 

(5 (Tr L/t [/„ - 1) = — y cff„„ exp (iTr T„„ - iTr C/t(T„„ ® 1jv)[/„) (B.IO) 

for n/m and n = m, as above. Just as the trace Tr , the direct product is only needed in the case (3 = 4. 

We order the columns U n , n = 1, . . . , N of the matrix U in a vector U = (Ui, U2, ■ ■ . , Un) t with iV 2 elements. 
For /3 = 1, 2, the elements are scalars, for ft = 4, they are quaternions. Collecting everything, we can re-write the 
measure (B.6) in the form 

d\U] f / \ 

= (27r) „> ( „l 1)/2 / d[T]exp(iTrT-iTr[/t (T8liv) ^ . (B .n) 

To use this in the integral (3.15) for the matrix Bessel functions (x, k), we also take advantage of the relation 

Tr U~ x xUk = TxU ] {x®k)U (B.12) 

which allows us to write 

r(/3) 



(iTiW(x®k-T® 1 N )U 



exp 1 « . 



(27T) 



AT 



d[T]exp(iTrT) 

Det~' 3/2 f a; ® - T ® l^v) • (B.13) 



Thus, the integration over U could be done as a Gaussian one and gave the result (B.l). Obviously, the Gaussian 
integrals over U only converge, if the diagonal elements of T have a proper imaginary increment. 

Formula (B.13) yields immediately the integral equation (B.3). Upon making the change of variables 

T = x l ' 2 T'x l l 2 , implying d[T] = Det i+^-D^ d[r] > (B14) 
we bring x into the exponential function and remove it from the determinant. We diagonalize T' = V'~ l t'V' and find 
Det^ 2 (x®k-T® ljv) = net~ pN/2 x _/S/a . (B.15) 

The int egra l over V is then just the integral definition ( 3 . 1 5| ) of the matrix Bessel function <3>^(x, t') and we arrive 
at Eq. ( §£f ). 
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The normalization constants remain to be derived. -Conveniently, they nicely relate to a special form of Selberg's 
integral which is given in Eq. (17.5.2) of Mehta's bool£3, 



j N = f d[t]\A N (t)\ 2 ^ f[ { ai + it n y bi (a 2 -it n y 

71=1 

(2n) N 



(a a +a 2 )( bl + 62 ) Ar -^ Ar ( Ar - 1 )- JV 

X 1 T(l + (n + l)7)r(6i + 62 - (JV + n - 1)7 - 1) 



n 



(B.16) 



r(i + 7 )r(6i - n 7 )r(6 2 - n 7 ) 

We now put x = or ft = and have ®$(0,k) = 1 or $^(x,0) = 1 on the left hand side of Eq. flB.13| ). We 

diagonalize T = V~ X W and use the invariance of the integral. Employing the measure ( |3.7| ) and the constant 
given in Eq. (3.E), we find the condition 



M {0) C {P) 0N/2 
_ N N n 



(2tt 



(B.17) 



We map this onto Selberg's integral (B.lq) by setting 7 = (3/2, h\ = [3N/2 and a 2 = 6 2 , by using 



lim - — — 

a 2 ^oo (a 2 - it n ) a2 



— exp(it ri 



(B.18) 



and by considering a^ 12 Jn in the limits a% — > and a 2 — > 00. With the help of some standard asymptotic formulae 
for the r function, we obtain M^P and, eventually, the constants and in Eqs. (|B.2|) and (B.4). 



APPENDIX C: FOURIER— BESSEL ANALYSIS 



The Fourier-Bessel Analysis involving matrix Bessel functions was discussed by Harish-ChandraE] in a general and 
formal way. To show the connection to our results, we summarize here some essential features of the Fourier-Bessel 
analysis on an explicit level. 

We write the Fourier transform of a function f{H) as 

F(K) = Df J d[H] exp (iTr HK) f(H) (C.l) 

where the matrices H and K have the same symmetries. If we choose a symmetric normalization, 

n(/3) _ I (no) 

N (2^)^/2^/3^(^-1)74 ' ^■ A > 

we can write the inverse transform as 

f(H) = D ( P ] J d[K]exp(-iTrKH)F(K) . (C.3) 

We notice that, according to Eq. ( 3.11 ), the Fourier transform of the constant Djp is the 5 distribution 5(K) and 
vice versa. 

If / is an invariant function such that f(H) = f(x), its Fourier transform turns out to be invariant as well, 
F(K) = F(k). Introducing eigenvalue-angle coordinates, we easily find 

F(k) = Dfc ( f J d[x]\A N (x)f^\x,k)f(x) (C.4) 

for the Fourier transform and 
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f(x) = D^cf J d[k]\A N (k)f$f*(k,x)F(k) 
for its inverse. We now insert the transform ( |C.4| ) into the inverse ( |C.5[ ) and conclude that 

(Dfc^Y f d[k]\A N (k)f$f(x,k)*f*(k,y) 



(C.5) 



det[S(x n - y m )] r 



-1,...,N 



\A N (x)A N (y)\0/* 

This is the analogue of Hankel's expansion of the 5 distribution. From Eq. ( |C.6| ) , the formula 

1 det[5(x n - y m )] n ,m=l,...,N 



dfi{U) S (?7 1 xU - y) = 



C 



09) 

AT 



\A N (x)A N (y)\^ 



obtains. To see this, we introduce a matrix G having the same symmetries as H and write 

5{H-G) = (D^) 2 J d[K] exp (-iTrK(H - G)) . 
Averaging over the diagonalizing matrix U of H yields 

J dp,(U)5(H-G) = (D^f J d[K\$yp*(x,k) exp(iTri^G) , 



(C.6) 



(C.7) 



(C.8) 



(C.9) 



by using the invariance of the measure. We now introduce eigenvalue-angle coordinates for K and do the integral 
over V, the diagonalizing matrix of K, 



Jdv(U)S(H-G) = (D^) 2 C« 3) Jd{k}\A N (k)f^*(x,k)$f(k,y), 



(CIO) 



where we have, once more, employed the invariance of the measure. Since the right hand side of this equ ation does 
only depend o n the eigenvalues y of G, we may replace G on the left hand side with y. Together with Eq. (C.6), this 
gives formula ( |C.7D . 

For the convolution in matrix space of two functions fi(H) and /2(H), we straightforwardly find the generalization 
of the standard convolution theorem, 



f(H) = J d[G]h{G)f 2 {H - G) = J d[K] C xp(-tTiHK)F 1 (K)F 2 (K) 



(C.ll) 



where G has the same symmetries as H. The functions Fi(K) and F-AK) are the Fourier transforms of fi(H) and 
/2(H), respectively. If the functions are invariant, the second of Eqs. (C.ll) acquires the form 



f(x) = Cf / d[k]\A N (k)f^ ) *(x,k)F 1 (k)F 2 (k) . 



(C.12) 



On the other hand, we find from the first of Eqs. (C.ll) 



f(x) = Cf J d[y] | A N (y)f h(y) f 2 (x,y) (C.13) 

where y are the eigenvalues of G. This formula is a convolution in the curved space of the eigenvalues. The second 
function is given by 



f a (x,y) = / dn(U)f 2 (x-U^yU) . 



(C.14) 



We insert the Fourier integral for Fi(k) according to Eq. ( |C.4j ) into Eq. ( C.12| ), compare with Eq. (C.13) and obtain 
the Fourier decomposition 



f 2 (x,y) = D^C^ ) ld[k]\A N (k)f^ ) *(x,k)F2(k)^\k,y). 
Formulae (|cj) and (|C/7|) can be viewed as special cases of these results. 



(C.15) 
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APPENDIX D: ACTION OF THE LAPLACEAN ON THE RADIAL FUNCTIONS FOR ARBITRARY 



We make the notation more compact by defining 



d/j,(x', x) — d/J,(x', x) exp I i I x n — x' n ) ) , (D.l) 

V \n=l n=l / / 



where the measur e is given in Eq. fl5.6p . To prove the identity (5.9), we write the integral using functions. The left 
hand side of Eq. ( |5.9| ) reads 

A, /^(x^a ; )^ 1 (x'J)[]e(x ^ -4)^ e ^■- a; ^H2:'] , (D.2) 

J i>j j>i 

where now the integration domain is the real axis for all variables. Thus, we can directly calculate the action of the 
operator A x onto the integral. We find 



A x f Jl(x' ,x)^l 1 (x' ,k)l[Q(x t ~ x'^llQi^ - Xl )d[x'} = 
J i>j j>i 

J ^(x'.fejjje^-^jje^.-xo 

i>j j>l 

( A i' ) +/> E TZ, \ V 2 - H(x\x)d[x'] 
\ n7 t m 1 ™ m) j 

+ J ^l 1 {x\^Ji{x\x)A x l[Q(x i -x' j )l[e(x' j -x ) d[x'] 

i>j j>l 

N 

+2 J ^Jj^^^ij^ne^-^nec^-^dM , (d.3) 



i>j j>l 

where we define the operator 



A ( 7 )= Y**- y p (J. 9 -\ . (d.4) 



dx^ x n x m V dx n dx m 

n—l n<m x 



By a series of integrations by parts, the operator AL acting on /I(a:', a;) is transformed to A^/ acting only on fc). 
At taking the derivative of the functions, we notice that only adjacent levels contribute, because otherwise terms 
like Q(xi — Xj) with i < j arise which annihilate the integral due to the chosen ordering. Therefore, we can write 

^-l[Q(x i -x' j )l[Q( X ' j -x l ) = 

i>3 j>l 

J|(©^nn' ,jHn-l)'n)(Hx n - x' n )&{x' n _ x ~ I ") 

- - x n )6(x n -x£j) , (D.5) 

where n(®^nn'> ®^(n-i)'n) is short-hand for the product on the left hand side of Eq. ( D.5 ) without the two factors 
Q{x' n _i — x n )Q(x n — x' n ). Importantly, this product is symmetric in x' n _ 1 and x' n . The second derivatives yield 

Q( Xi - X' 3 ) Q{ X ' 3 - Xl) = Y[(e^nn>, O^n-iyn) 
OXn i>3 j>l 

(5'(x n - x' n )Q{x' n _ 1 - x n ) + 5' {x' n _ x - x n )Q(x n ~ x' n ) + 

- x n )S(x n - <)) . (D.6) 

The last term vanishes upon integration, since it is symmetric in x' n _ 1 and x' n , whereas the rest of the integrand is 
antisymmetric due to the Vandermonde determinant An-i(x') in the measure (|5.6|). Differentiation with respect to 
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— J{Q{x i -x' j )J{@{x' j -x l ) 



i>3 



J>1 



IJ(0^n'(n+l)) @^nn>)(S(x' n - X n+ i)Q(x n - x' n ) 

- S(x n - x' n )Q(x' n - x n+ i)J . 



(D.7) 



Integration by parts of the first term of the right hand side of Eq. (D.3) yields 
A X J ftx^x^^ix'^lleixi-x^neix'j-xOdix'} = 

i>j j>l 

J Ji{x\ x)A x ,^l 1 (x', k)d[x'\ -k 2 N J Jl{x\ x^l^x', k)d[x'] 
\S(x n - x n )Q(x' n - x n+ x) + S(x' n - x n+ i)Q{x n - x' n )\ 



d d 1 v-^ p 1 ^ P \\ , ... ,. 

OXn OXy, L Xn Xm , X^ X m J J 

>■■■■-■■■■ m^n n " L J J 



(D. 



Inserting in Eq. (D.8) the function /j,(x',x) as given in Eq. (DA) and Eq. (5.6) we find after a straightforward 
calculation 

A X J Ji(x\x)*% ) _ l (x',k)d[x'] = 

Jl(x',x)A x ,^ ) _ 1 (x',k)d[x'] - k 2 N j ]i{x' ^x^l^x' ,k)d[x'] 



with 



+2 J ^l^X^k)^ ^n{®?n>( n +l),®*nn>)(^(x' n ]y n ;X,x') 

(s(x n - x' n )Q(x' n - x n+ i) + 5{x' n - x n+ i)<d{x n - x' n ) jj Jl(x', x)d[x'] 

( N - 1 I i \ 

g(x n ;x' n ;x, x') = (P/2 - 1) £ — 3" - E — ~ 

i a, n x m ^ n •km f 

5 «;x„ ;a ;,,') = (/3/2-l)(E^r-E^^] ■ 

V m^n 11 "* m=l '* / 



(D.9) 



(D.IO) 



We now can perform the integration of the 5 distributions in Eq. ( DS, ) . We notice that the difference (g (x' n ; x n ; x , x' ) — 
g(x' n ;x n ;x,x')) vanishe s lin early, whenever x' n approaches one of the boundaries of its integration domain. Thus the 
second integral in Eq. (D.9) yields zero as long as the measure diverges less than (x n — x^) -1 when x' n approaches 
x n - This is always the case for p > 0. Collecting everything, we arrive at the identity (|5.9|). 



APPENDIX E: SYMMETRY OF THE RADIAL FUNCTIONS FOR ARBITRARY 



Applying the recursion formula (5.5) to all N — 1 levels, we can extend Eq. flf.5|) to arbitrary p and write 



. N-l 

$f ( x ,k)= j n ^o* (n) >s (n-1 
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N—n 

E 

rn — 1 



» u 

,™ I ft, 



N— n+1 



(E.l) 



where x*- - 1 = a;. Analogously, we also find 



. JV-l 

< } (fc^)= / n dn{k^M n ~ x 

J n=l 



exp 



exp 




N — 7 



'0 



%N-n+l 



m— 1 



(E.2) 



with = fc for the solution of the different ial e quatio n w hich results from Eq. ( |5.f|) by interchanging x and fc. We 
have to show that these two radial functions (E.l) and (E.2) are identical. To this end, we change in Eq. (E.l) on the 
n th level the variables x 



(n) 



,(N-n) to kt\ m 



N-n-1, (n-1) 



1 = 1 i X 



1, . . . , (N — n) by setting 



(n-l) 



(«-l) 



-fc^) 



(E.3) 



for ri = 1, ... , (JV — 1). These are, on the n th level, N — n + 1 equations for making a change of N — n variables. 
However, one has 



N-n+i 

E ' 

m— 1 



„(«) 



(E.4) 



on all levels which eliminates o ne o f the N — n + 1 equations. 

Of course, the substitution (EJ5) is motivated by the radial Gelfand-Tzetlin coordinates which we introduced to 

(n) 

construct the recursion formula for f3 = 1, 2, 4. In this case, the r m are the moduli squared of a c olum n of a matrix 
U G U(N — n; (3). Here, we do not use this connection to ma trice s and groups. We simply view Eq. (|E.3|) as a standard 
change of variables in an integral. We underline that Eq. (EJ3) does not involve (3 at all. Equation (1^4) is just the 
normalization o f a column of U for (3 — 1,2,4. Since it is independent of (3, it also holds for arbitrary (3. One can 
also verify Eq. (E.4) by a direct calculation. 

The original domains of integration are Xm ^ < x*m < x^, \ . In these boundaries, the are positive definite. 

Hence, to satisfy this when changing the variables, we must have km ^ < fcm'* < ^m+i ^ or ^ ne new domains of 
integration. 



To work out the measure in the new variables km , we interpret Eq. (E.3) as a change to the integration variables 



Tro , too. This yields immediately 

A N _ n (x^) 



A 



N-n+l 



(x( n - 



■ d[x 



(n)l 



dix{r {n) ) 



(E.5) 



The first equality sign goes back to the radial Gelfand-Tzetlin coordinates. We may use this piece of information, 
because it is independent of (3. The second equality sign is simply due to Eq. (E.3). Using this result, we find for the 
full and j3 dependent measure 



dfi(x 



(«) X (, n -V^ 



G 



(P) 

N-n+l 



'iL 



>-l) 



W-2)/2 



A2 ( r (n-l)\ 



(a,(n-l)) 



±N-n+l 
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(«)l 
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(/3-2)/2 
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(n-l) 



*f B) ) 



A 2 
N- 



.n+i(fc (n - 1) ) ; 



Aiv-n+iCA;^- 1 )) 



(/3-2)/2 



dn{k {n \k {n -^ 



(E.6) 



For /3 = 1, 2, 4, this result is a direct consequence of the invariance of the group measure dfx(U). Here, we have derived 
it for arbitrary [3. This, in turn, implies that the invariance of the group measure dfi(U) is embedded into and reflects 
much more general features. 

We now collect all these intermediate results and plug them into Eq. (E.l). Apart from the expressions in the 
exponential functions, we ha ve fu ll agreement with the right hand side of Eq. (E.2). Hence, it remains to be shown 
that the change of variables leads to the identity 




„(JV-i) 



N-r 



4 n) ZiV-n+l + k[ N - 1] Xl 



(E.7) 



m— 1 



Since th e sy mmetry relation (5.3) holds for /3 = 1, 2,4, we know that Eq. (E.7) must be true in these cases. However, 
as Eq. ( |E.7| ) doe s not involve (3 a t all , it must also be valid for arbitrary (3. Inserting this into the right hand 
side of Eq. (E.l), we recover Eq. flE.2| ), as desired. We notice that this line of arguing cannot be spoiled by any 
other contribution to the argument of the exponential functio ns, b ecause all other terms in the integrand are purely 
algebraic. This completes the proof of the symmetry relation (5.3) for arbitrary (3. 



APPENDIX F: CALCULATION OF THE NORMALIZATION CONSTANT 

In the previous App. [e], we introduced the coordinates r' n = r„ , n = 1, . . . , N on the first level of the recursion. 
They are the moduli squared of the coordinates on the unit sphere in the complex N dimensional space. Thus, it is 
natural to use the following type of hyper spherical coordinates 

n-l 

' n = cos d n Y\_ sm &u , n = 1, . . . , (N - 1) , 

r N = sin#jv_i sintf,, (F.l) 
i/=i 

where the positive semidefiniteness of the r' n restricts the domain of integration to < $ n < 7r/2, n = 1, . . . , (N — 1). 
Thus, we integrate over a (2 N ) th segment of the unit sphere. The measure 

N-l 

Mr') = ]J sin 2(JV " n)_1 i?„ cos0„d0„ (F.2) 

n=l 

is, apart from the phase angles, the measure on the unit sphere. Collecting everything, we have 

(N \ 

N-l r7T /2 

= G { V TT / sin^-)' 3 - 1 K cos' 3 ' 1 d n dK 

n = l J 

(0) "jj 1 T((N-n)(3/2)T((3/2) m T N ((3/2) 

N JLJL 2T((N -n+l)(3/2) N 2 N - 1 T{N(3/2) [ 1 



where the integral over d n is just Euler's integral of the first kind. 
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APPENDIX G: TRANSLATION INVARIANCE OF W$l(X,K) 



We shift every x n in the the recursion formula ( |5.5| ) for arbitrary (3 by a constant x and obtain 
®jP(x + x, k) = y dfi(x',x + x) exp x + iVx - ^ sc'J few^ 



(G.l) 



with x n + x < x' n < x n -\-i + x as the domains of integration. The change of variables x' n — > x' n + x removes x from 
the measure given in Eq. (El]) and the domains of integration, we find 



N N-l 



^l^x' + x,k) . 

We want to employ an induction. We assume that the radial functions for arbitrary have the property 



^\x + x,k) = exp Hx^k^j ^\x,k) 



(G.2) 



(G.3) 



If this is correct for N — 1, formula (G.2) implies that it is also true for N. T he in duction starts with N = 2 where 
the correctness of Eq. ( G.3 ) is immediately obvious from the explicit solution ( 3.21 ) for arbitrary f3. Thus, Eq. ( |G.3| ) 
is valid for all N . 

Since the k n are arbitrary and since the sum over all k n is invariant under the permutations oj(k), the property ( |G.3| ) 
must also be true for every function - u {,x, k) with oj G Sn- We compare this with the expression 



N 



®N,u>( X + x i k ) = ex P ix fc " 



exp 



o(n) 



\A N (x)A N (k)\W 



W { N l(x + x,k) 



(G.4) 



which results from the Hankel ansatz (5.1C). Hence, we conclude that we necessarily have 

Wj£l(x + x,k) = Wghx,k) . 



(G.5) 



This is the translation invariance. 



APPENDIX H: CALCULATION OF ^ ] {X,K) 



W e per form the calculation for §^ (—ix, k) to avoid inconvenient factors of i. The operator L x 



Eq. (5.13) splits into two parts. The first part 



defined in 



N & 

Y — 



(H.l) 



does not change the order in k, while the second one, 



A 



N 

2 Vfc„ 



1 dx n 



(H.2) 



raises the order in k by one. Since we can restrict ourselves to one element of the permutation group, we discuss 



only the identity permutation in the sequel. The symmetry of x and k together with the result for <i> 3 (x, k) suggests 
one to try an expansion in the composite variable Zij as defined in Eq. (6.i 
symmetric functions 



To this end we define the elementary 
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= e n 

»1J1<*333<— <ti/jw 1 = 1 



■H3i 



(H.3) 



Here, we assume the following ordering of the composite index {iiji}, ii < ji- We say {iiji} < {i m jm} if h < *m or 
= i m and ji < j m . All indices run to N. The highest o rder elementary symmetric function is of order N(N — l)/2 
and is given by A]v(a;)A7v(fc). The asymptotic formula ( |5.14 ) yields the leading term for large arguments. It is the 
starting point for a recursion in powers of z . 



N(N-l)/2 



u=0 



(H.4) 



where y y (z) i s a s ymmetric function of order v in x% and fej. We investigate the action of the two operators defined 
inEq. (|HT|) ,B) and find 



A^ fc e„(z x ) 
A a; ,fce ly (z _1 ) 
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2^ 
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2y-l(^ Ilm ) 



(H.5) 



e "( Z inra) 
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^ n k 
^ m k 



(H.6) 



The function e„(z^„ m ) is the elementary symmetric function e„(z) with all terms con tain ing z nm omitted. For 
^ = 0,1,2 we simply have p u (z^ 1 ) = (— 2) ly e„(z~ 1 ). For i/ > 3 the last term in Eq. (EL5) causes corrections to 
the elementary symmetric functions. This arises due to the mixed derivatives which have to be taken into account 
in the action of A x ,k onto e l/ (z _1 ) for v > 3. Because of this term the Hankel Ansatz seems becomes increasingly 
cumbersome as higher values of N are considered, since more and more correction terms have to be constructed. So 
far, the construction was only possible for N = 4. To construct the correction terms explicitly for the case N = 4, we 
define a new set of symmetric functions as follows 



N 

Z kl Z km Z lm e v-3( z ) 



/ j "kl ^km^lm 
k<l<m 



(H.7) 



Again we have to investigate the action of A^ fe and A x ^ k on /^(z *). We find 

1 



N 

E - 



2 fs^^nm) 



and 



N 



Ax,kh{z *) 



o V - z- 1 z- 1 

(~ _ t XI nk z mk i 



(H. 



(H.9) 



thus /^(z 1 ) is the desired correction term. We have 



(H.10) 



Fortunately, due to Eq. (H.8) in the next step the correction term itself has not to be corrected and we find 

Piiz- 1 ) = 2 4 (e 4 (z- 1 ) + ij- 4 (z- 1 )) . (H.ll) 
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Up to now the se re sults are valid for arbitrary N. The action of A Xt k onto the symmetric function f±[z 1 ) is not as 
simple as Eq. ( H.8 ). After a series of manipulations we arrive at 



JY 



A x , k hiz- 1 ) = -4 Y, 



7Z _ T \2 f^ Z ^nm) 



N 



2 Z nk Z mk fc( Z ) 



k^r. 



j^mk 



(H.12) 
(H.13) 



The contribution ( |H.6| ) has to be added to this expression stemming from the action of A x .k onto 64(2 1 ). On the 
other hand we calculate 



iV 



Ax, k f 5 (z x ) = -2 ^ - _ X / 4 (z^J 

l^n x in J 



AT 
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, ,. ^2 Z nk Z mk e 2(z/ nm ) 
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(H.14) 
(H.15) 



Thus, we have to find yet another correction term to compensate the second term in Eq. (H.12 ). We define 



E IT Z iyij E Ziri i ~ E Z jl Z jm Z kl Z km Z lm 



(H.16) 



i<k 

Km 



and see that A^^/^z -1 ) yields exactly the desired second term of Eq. (H.13). Pushing forward this procedure 
becomes more complicated step by step. There seems to be no obvious way of constructing the additional terms. 
Apparently for higher orders the correction terms also involve an increasing amount of indices. Nevertheless for N = 4 
we are already at the end of the recursion. Then the general expression 



P5 



reduces to 



p 5 (z- 1 ) = -72e 5 (z- 1 ) 



(H.17) 



(H.18) 



The last step can readily be done, since the action of A x ,k onto es(z x ) is already known by Eq. (H.6). Thus we 
arrive at 



p 6 (z- 1 ) = 288e 6 (z~ 1 ) 



Importantly, we have 



A^fe e 6 (z 1 ) = A X}k 



1 



A 4 (a;)A4(fc) 







(H.19) 



(H.20) 



That means, the sequence finishes after the sixth step. Collecting everything and observing that, for N = 4, /s(z) = 
2e§(z) and fe(z) = 4eg(z), we obtain 



Wi 4) (x, k) = ^(-2)^e,(z- 1 ) + ^(-2)"- 1 / t/ (z- 1 ) - 8 e 5 (z- 1 ) + 96 e^ 1 ) 



(H.21) 



v— 1 v— 3 

This can be rewritten more in a more compact way as 



af'(x,k) 
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5 e n (2-««)+zE n p-^) 



(H.22) 



l<m<n i <i 



Kir 
fe<r. 



i<J ^ifc 



which yields Eq. (6.1C) 
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